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Abstract 

We study the BGK approximation to first-order scalar conservation 
laws with a flux which is discontinuous in the space variable. We show 
that the Cauchy Problem for the BGK approximation is well-posed and 
that, as the relaxation parameter tends to 0, it converges to the (entropy) 
solution of the limit problem. 
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1 Introduction 

In this paper we consider the equation 

dtf + d,{k{x)a{On^^^^^^^, t>0,xeR,^eR, (1) 
with the initial condition 

rit^o -/o, inM, xMj. (2) 

Here k is given by 

k = fcLl(_oo,0) + fci?I(0, + oo), 

where 1b is the characteristic function of a set B, ^ M- a(^) is a continuous 
function on R such that 



Vw e [0,1], [ a{^)dC > 0, [a 
Jo Jo 



iOd^ = 0, (3) 
and, in (IT|), Xw , the so-called equilibrium function associated to /' is defined 

by 

U^t,x)^ I rit,X,Odi, Xa(e) = I]0,a[(O-I]a,0[(O, 
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for t > 0, a; e R, f e R, a e M. 



Eq. ([Tl is the so-called BGK approximation to the scalar conservation law 

dtU + d^ikix)A{u))=0, A{u)= f aiOdC (4) 

Jo 

The flux {x,u) n- k{x)A{u) is discontinuous with respect to a; G M, actually (H]) 
is a prototype of scalar (first-order) conservation law with discontinuous flux 
function. In the last ten years, scalar conservation laws with discontinuous flux 
function have been extensively studied. We refer to the paper ^BKQSj for a 
comprehensive introduction to the subject and a complete list of references. Let 
us simply mention that the discontinuous character of the flux function gives 
rise to a multiplicity of weak solutions, even if traditional entropy conditions are 
imposed in the spatial domain apart from the discontinuity. An additional cri- 
terion has therefore to be given in order to select solutions in a unique way. For 
scalar conservation law under the general form dtu + dx{B(x,u)) = 0, where 
the function B is discontinuous with respect to x, several criteria are possi- 
ble |AMG05| . For B{x,u) = k{x)A{u) as above, the choice of entropy solution 
is unambiguous (see [AMGOSj . Remark 4.4) and we consider here the criterion of 
selection first given in |Tow01j . A kinetic formulation (in the spirit of LPT94j) 
equivalent to the entropy formulation in |Tow01] has been given in [BV06j . In 
particular, solutions given by this criterion are limits (a.e. and in L^) of the 
solutions obtained by monotone regularization of the coefficient k in ([¥]), e.g. 

ke{x) = kLlx<-e{x) + I — X H 1 l^^<x<e + kRle<x, £ > 0. 

The kinetic formulation of scalar conservation laws is well adapted to the analy- 
sis of the (Perthame-Tadmor) BGK approximation of scalar conservation laws. 
Developed in jPT91) . this equation is a continuous version of the Transport- 
Collapse method of Brenier |Bre81[ IBre83| . BGK models have also been used 
for gas dynamics and the construction of numerical schemes. See for example 
the book of Perthame [Per02j for a survey of this field. 

Our purpose here is to apply the kinetic formulation of [BV06) to show the 
convergence of the BGK approximation. To this aim, we first study the BGK 
equation in itself in Section [2j In Section [3l we introduce the kinetic formula- 
tion for the limit problem. We also introduce a notion of generalized (kinetic) 
solution. Definition ini We show that any generalized solution reduces to a mere 
solution, i. e. a solution in the sense of Def. S) This theorem of "reduction" is 
Theorem [71 Then in Section HJ we show that the BGK model converges to a 
generalized solution of (|4]) and, using Theorem[71 deduce the strong convergence 
of the BGK model to a solution of (g]). Theorem [TTl 

A key step of the whole proof of convergence is the result of reduction of The- 
orem [T] Its proof, given in Section 13. 2[ is close to the proof of uniqueness of 
solutions given in |BV06| . A minor difference is that we deal here with gene- 
ralized solutions instead of "kinetic process solutions". There is also a minor 
error in the proof given in |BV06j (specifically, the remainder terms Ra,e,s and 
Qp,v,a in Eq. (PT|) and ([5^ of the present paper are missing in |BV06] ). We 
have therefore given a complete proof of Theorem [71 
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We end this introduction with two remarks: 

• the BGK model provides an approximation of the entropy solutions to ^ 
by relaxation of the kinetic equation corresponding to Q. A relaxation 
scheme of the Jin and Xin type applied directly to the original equation 
dH) has been developed in |KKR04) . 

• in the last chapter of |Bac05| . is derived the kinetic formulation of scalar 
conservation laws with discontinuous spatial dependence of the form dtu + 
dx{B{x,u)) = (which are more general than ([4])). We indicate (this 
would have to be proved rigorously), that in case where our approach 
via the BGK approximation was applied to this problem, the solutions 
obtained would be the type of entropy solutions considered in |KRT03| . 

Notation For p,q € [1, +oo], we denote by L^.L'^ the space LP{M.a:; L''(K.j)) and 
by LjlP the space {R^; Lp (R^)) . 

We also set sgn+(s) = I{s>o}, sgn_(s) = -l{s<o}, sgn = sgn+ + sgn_, s G M. 



2 The BGK equation 

2.1 The balance equation 

By the change of variables /^(t,a:,0 — ei f^{t,x,^), Eq. ([T]) rewrites as the 
balance equation 

£ 

t 

with (unknown dependent) source term ■ Hence, we first consider the 

following Cauchy Problem for the balance equation: 

dtf + dx{k[x)a{Of)^g. t>0,xeR,eeM, (5) 

/|t=o =/o in X M^. (6) 

Proposition 1 Suppose that kn ■ > 0. Then Problem (^-^ is well posed 
in L\LP., l<p< +oo; for all fo £ L^LP, T > and g £ L^{]0,T[; LILP), there 
exists a unique f £ C{[0,T]; L^LP) solving ^ in V {]0,T[xRx x R^) such that 
/(O) = /o- Besides, we have 

< Mk f II/oIIljls + r hm^iL^js] , (7) 



where Mk =max(|^,|^). 



Proof: Since (O is linear, it is sufficient to solve the case 5 = 0. The general 
case will follow from Duhamel's Formula. Assume without loss of generality 
kR, kL > 0. Let A+ := £ R; a(C) > 0}. Then, for fixed £, & A+, and although 
/c is a discontinuous function, the O.D.E. 

X{t,s,x,^) = k{X{t,s,x,0)a{^), teR, (8) 



3 



with datum X(s,s,x,^) = x has an obvious solution for a; ^ 0, given by 
s, X, ^) = X + — s)kfia{S,), t > s, when x > 0, and by 

X + {t - s)kLa{S,) if t<s 



when X < 0. Denoting by s+ = max(s, 0), = s+ — s the positive and negative 
parts of s e R, and introducing 



afc(x) = I{x>o} + ^I{x<o}, 



this can be summed up as 

X{t,s,x,(,)^{ak{x)x + {t-s)kRa(Cj}+ -{x + [t-s)kLa{Cjy, t > s. (9) 
Similarly, we have, for the resolution of ^ backward in time, 

X{t,s,x,^) = {x + {t-s)kRa{()}+ -{Pk{x)x + {t-s)kLa{()}-, t < s, (10) 
where 

Pk{x) = -r—^x>o} + I{.T<o}- 

A similar computation in the case a(^) < gives the solution to ([S]) by ([S]) for 
{t - s)a(C) > 0, ^ for (t - s)a(^) < 0. For the transport equation {dt + 
k{x)a(£_)dx)f* = 0, interpreted as 

this yields the solution 

^*{t,x,0=i^{X{T,t,x,0,0. 

which satisfies the terminal condition >p*(T) = ip. We suppose in what follows 
that ip is independent on compactly supported and Lipschitz continuous. 
Then, a simple change of variable shows that, for every t G [0,T], for a.e. ^ G M, 

y*{t,;0\\Ll<Mk\mLi, Mfe=max(^^,^^,l<9<+oo. (11) 

If / e C([0,T]; L|LP) solves dS])-®, then, by duality (note that (p* is Lipschitz 
continuous and compactly supported in x if "0 is) we have, for t G [0, T], for a.e. 

/ fiT,x,Oi^{x,Odx^ f /o(x,O^*(0,x,Orf2:. (12) 

In particular, the estimate (jlip where q = conjugate exponent of p gives, for 
a.e. C e M, 

||/(T,-,OIIlS <^'^fcll/o(-,OllL£, 
and then by Duhamel's principle, for g 0, 

||/(T,-,OIL£ <Mfc (||/o(-,OIIl£+ rMt,;mL^Jt] . (13) 
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The estimate (O and uniqueness of the solution to ([S])-® readily follows. Exis- 
tence follows from (pi)- (|10l) - p^ . from which one derives the explicit formula 

the coefhcient J{t,x,^) being given by 

kL 

J{t,X,^) = I{x<0}U{a;>tfc„a(C)} + 7~ I{0<2;<tfeHa(J) } 

if a(^) > and 

J{t,X,^) = l{x<tkLa{^)}u{x>a} + ■^I{tfeta(e)<2;<0} 

if a(0 < 0. 

2.2 The BGK equation 

Denote by T(t)/o the solution to ([S])-® with g = 0, i.e. 

r(i)/o(x,o = j{t,x,Ofo{x{o,t,x,o,0, 

X given by ®-(|ini)- 

Definition 2 Let fo e ^^(R:^ x R^), T > 0. yl function G C{[0,T]; L\Rx x 
R^)) is said to be a solution to (^-(^ if 

nt) - e-*r(t)/o + - e-iT{s)xu^f^t-s)ds, ^ f fi^d^, (14) 
for all t e [0,T]. 

Theorem 3 ^sstime fcj^ • fc^ > 0. Lef /o G i^ll^a: x R^), T > 0. There exists 
a unique solution f^ G C{[0,T]; L^(Rx x R^)) to Denoting by Se{t)fa 

this solution, we have: 

1. \\{SS)fl " SS)flr\\L^WL.yM,) < Mfe||(/5 - /o')+||li(K.xE,) 

2. 0< sgn(0/o(a;,0 < 1 a.e. ^ < sgn(05,(t)/o(x, < 1 a.e. 

5. ?//o - Xu„, uo G L°°(M), < uo < 1 a.e. then < 5e(i)/o < Xi- 

Proof: the change of variable {t' , x') = e{t, x) reduces ([T]) to the same equation 
with e — I. We then have to solve / = F{f) for 

F{f){t) e-*r(t)/o + e-'T{s)xuit-s)ds, u = fi^d^. 
By d?]) and the identity 

J \Xu-Xv\{£,)d^=\u--v\, u,i;GR, 
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we have F: C([0,T];ii^) C([0,r];Li^) and is a (1 - e"'^) contraction 
for the norm 

11/11 = sup ||/(i)llLi(R.xRe)- 
te[o,T] 

Indeed, we compute, 

im^rn ^ Fif'mhl^ <J^e-^\\T{s){XuHt^s)~Xunt-s))\\LlJs 

= / (^^''WXuiit-s) - Xu'>{t-s)\\L^ ds 

Jo 

s) — — s)\\]^ids 

By the Banach fixed point theorem, we obtain existence and uniqueness of the 
solution to (H])-©. Since < sgn{^)xu{0 ^ 1 ^-S- we have 

0<sgniOF{f){t,x,0 < 1 a.e. 

if < sgn(f )/o(a;, ^) < 1 a.e. This proves the point 2. of the Theorem. The 
point 1. follows from the inequality 

^sgn+(/-g)(g(/)-g(.g))de<0, f,geL'iR^), Q{f):^Xjfde~f, 

that is easy to check, and from the identity 

f{t) = r{t)fo+ I T{s)Q{f){t - s)ds 
Jo 

for the solution to ([T])-®. If /o ~ Xua, < wq < 1 a.e. then = xo < /o < Xi- 
Hence the item 3. follows from 1. and the fact that any constant equilibrium 
function Xa, a G R is solution to ([1]). ■ 

3 The limit problem 

Assume /o — Xuo with uq G L°°{M.), < uo < 1 a.e. Set 

Aiu)= r aiOlio.i]iOdC (15) 
Jo 

Note that by ([3]), we have A > and A vanishes outside the interval [0, 1]. We 
expect the solution to ((I])-© to converge to the solution u of the first-order 
scalar conservation law 

dtu + d^{k{x)A{u)) ^ 0, t>0,x&R, (16) 

with initial datum 

u(0,x) = ^0(2;), xeR. (17) 
For a fixed T > 0, set Q =]0,T[xM^. 



Jo 
Jo 

< f e-'ds\\f^^ 
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Definition 4 (Solution) Let uq E L°°{M.), < < 1 a.e. A function u 6 
L°°[Q) is said to he a (kinetic) solution to li6]j - |77P if there exists non-negative 
measures m± on [0, T] x R x R such that 

• m-(_ is supported in [0,T] x Rx] — oo, 1], m_ is supported in [0,T] x R x 
[0,+K)[, 

• for all e C^{[0,T[xR x R), 
h±{dttp + k{x)a{£_)dxtp)d£,dtdx 

+ [ [ ho,±i>{0,x,Od^dx-{kL-kR)^ [ [ a{OHt,0,Od^dt 

d^^jdm±{t,x,Cj (18) 
where h±{t,x,Cj = sgn±(u(i,a;) - Cl, ^o,±(a^,0 = sgnj^{uo{x) - £). 



'Q 



'Q 



Proposition 5 (Bound in L°°) Let uq e L°°(R), < ito < 1 a.e. If u E 

L°°[Q) is a kinetic solution to i_? ^) -( |T7| ), then < m < 1 a.e. 

Proof: Consider the kinetic formulation ([T5|) for with a test function 

= ip{t,x)^i{£). 
If pi is supported in ]1, +oo[, two terms cancel: 

JR JR JK JR 

and 

d^ipdm+{t, x,£,) — 



by the hypothesis on the support of rn+. Hence we have 
h+{dtip + k[x)a{X)dx^)ti(C)d£,dtdx 



- {kL - kR)+ f f a{Ov{t, 0)md^dt - 0. 
Jo Jm 

A step of approximation and regularization shows that we can take /i(^) = I{>i 
in this equation. Since 

/ + 00 
a(C)df = A{+oo) - A{1) =0-0 = 0, 

and 

h+it,x,^)d^ ^ / li<u{t,x)d^ ^ iu{t,x) - 
1 ^1 
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/ + 00 n + oa 

/u(t.x) 
a{Od^ = sgn+(u(i, x) - l){A{u{t, x) - A{1)), 

we obtain 

{u- l)+dtifi + fc(a;)sgn_^(M - l){A{u) - A{l))d^ipdtdx = 0. 

It is then classical to deduce that (m — 1)+ =0 a.e. (see the end of the proof of 
Proposition [TUl after ([M)) ). i.e. u < \ a.e. Similarly, we show u > a.e. ■ 

Our aim is to prove the uniqueness of the solution to P^ - dTT)) . Actually, more 
than mere uniqueness of the solution to ([T6)) -([T7 | . we will show a result of 
reduction/uniqueness (see Theorem [T]) of generalized kinetic solution. To this 
purpose, let us recall that a Young measure Q — R is a measurable mapping 
(<:, x) v-^ Ut,x from Q into the space of probability (Borel) measures on R. The 
mapping is measurable in the sense that for each Borel subset A of R, (<, x) n- 
vt,x{A) is measurable Q ^> R. Let us also introduce the following notation: if 
/ e L^{Q X R), we set 

This is consistent with the notations used in Def. |3]in the case f — Xu- 

Definition 6 (Generalized solution) Let uq £ L°°(M), < mq < 1 a.e. A 
function f G iy^(QxR^) is said to be a generalized (kinetic) solution to U6]) -[T1^ 

/ <i Xi a.e., — 9^/+ is a Young measure Q — > M, 
and if there exists non-negative measures m± on [0,T] x M x R .such that 

• m+ is .supported in [0,T] x Rx] — oo, 1], m_ is supported in [0,T] x R x 
[0,+oo[, 

• for all il) £ C^{[0,T[xR x R), 

/ / f±idti^ + kix)a{Odx^)d^dtdx 
Jq Jm 

+ 11 fo,±^(.0,x,Od^dx-{kL-kHf [ [ aiOHt,0,Od^dt 

JmJm Jq Jm 

= 11 d^'ipdm±{t,x,^) (19) 
Jq Jr 

where fo,±{x,^) = sgn±(uo(a;) - 0- 

Theorem 7 (Reduction, Uniqueness) Let uo G i°°(R), < uq < 1 a.e. 
Problem il6 ]) - [Tl^ admits at most one solution. Besides, any generalized so- 
lution is actually a solution: if f £ L^{Q x R^) is a generalized solution to 
(03)- (T^), then th ere exists u £ L^{C^^ such that f — \u- 

To prepare the proof of Theorem [71 we first have to analyze the formulation 
(|19p and the behavior of / at t = and x = Q. 
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3.1 Weak traces 

Introduce the cut-ofF function 

LUeis) = / Peir)dr, p,{s) = e-^pie-h), s £ M, (20) 
Jo 

where p S (M) is a non- negative function with total mass 1 compactly sup- 
ported in ]0, 1[. We have the following proposition. 

Proposition 8 (Weak traces) Let f e L°°{Q x M^) be a generalized solution 
to IM-^n^- There exists fl° G L'^(R x M), F± e L^{]0,T[xR) and a sequence 
(?7„) i such that, for all ip £ Ll{RxR), for all 9 e Ll{]0,T[xR) (the subscript 
c denotes compact support), 



f±{t,x,0^'r,St)v{x,Odidtdx / / f^^{x,Ov{x,Odidx, 

■/R JR 

(21) 

f±{t,x,0k{x)a{(,)uj'^Jx)9{t,(,)didtdx ^ [ [ F±{t,^)e{t,Od(dt (22) 



as n +CXD. Besides, there exists non-negative measures rrJ± , rn± on and 
[0,r] X R respectively such that: 

• (resp. fn^) is supported in Rx] — oo, 1] (resp. [0,T]x] — oo, 1]), vrJ° 
(resp. fh-) is supported mR x [0,+oo[ (resp. [0,T] x [0, -|-oo[j, 

. for all if e C^{R^), 9 G C^{[0,T[xR), 

f fi'^dxd^^f fo,±vdxd^~ f d^^dm^^{x,0, (23) 

JR2 JR2 Jjj2 

[ [ F±ed^dt ^-{kL- kn)^ [ [ a{£_)0d£.dt 
Jo Jr Jo Jm 

T I- 

d^edfh±{t,0- (24) 



Proof: The first part of the proposition does not use the fact that / is solution. 
Indeed, since |/±| < 2, we have 



f±{t,x,CW^{t)dt 



■T nT 

<2 I \u'^{t)\dt^2 p,j{t)dt<2, 
Jo 



for all {x, ^) G R^. This gives in particular a bound in L'^{K), K compact of 
on f±{t, ■)uj'^{t)dt, hence existence of a subsequence that converges weakly in 
[K] . Writing R^ as an increasing countable union of compact sets and using 
a diagonal process, we obtain (PT|) . The proof of (|22l) is similar. To obtain (1^51) . 
apply the formulation to ^{t,x,£^) = (p{x,^){l ~ujjj^{t)). We obtain (^5)) 
by using (PT|) and setting 



ipdni^ {x, 5) 



lim 



V{x,0i^ - ^vn.i't))d'm±{t,x,0 
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for all non- negative G Cc(K^): the limit is well defined since the argument is 
monotone in n and it defines a non- negative functional on Cc(ffi^) which is repre- 
sented by a non-negative Radon measure. Similarly, applying the formulation 
(HH) to ^(t,x,i) = Q{t,i){\ - w^„(a;)), we obtain ^ with 

/ fedfh±{t,0^ lim / -w^„(a;))dm±(t,a;,e) 
for ah non- negative 6 £ Cci[0,T] x R). ■ 

Remark: Since < / < xii (EU shows that , resp. is supported 

in Mx] — cxD,l], resp. R x [0, -|-cxd[. Similarly, resp. F_, is supported in 
[0, T] x] — cxo, 1], resp. [0, T] x [0, +oo[. We use this remark to show the following 

Corollary 9 For all ip^ e L°°(M.^) supported in [~R,R\ x [-i?,+(X)[ (R > 0) 
such that d^ip- < (in the sense of distributions), we have 



lim / / f+Lu' {t)Lp_{x,i)d£,dtdx> / fo +ip-dxd^. (25) 
"■'^+°°JqJr 7r2 ' 

For all 9- € L°°{]0,T[xR) supported in [0,T] x [-i?, +oo[ (R > 0) such that 
d(^9- < (in the sense of distributions), we have 

lim / / f+k{x)a{Ooj' {x)e^[t,i)d£,dtdx > -{kL-kR)+ [ f ai^O-d^dt. 
^-^"^ Jo Jr 

(26) 



n— >-t-oo 



Q JR 



Proof: Note first that each term in (|25p is well defined by the remark above 
and that, by (pij) . 



lim / f+{t,x,Oiolit)^_{x,Odidtdx= fl\^d^dx. 
By regularization (parameter e) and truncation (parameter M), we have 
if? - fo,+)^-dxd^ = / (/T" - fo.+W'l'" dxd^ + ,y(e. A/), 



where lim rile.M) = 0. More precisely, we set 

^-^ = iV- *A) X XM, 

where ip^ is a (smooth, compactly supported) approximation of the unit on R^ 
and xm is a smooth, non-increasing function such that xm = 1 on ] — oo, M], 
XM = on [M + 1, +oo[. Apply ^ to ip"^^^ to obtain 

(/? - /o.+)¥'^rf^de = - / 94/j^'dm;« (x, + v{£, M). 

JR2 

For A/ > i? + 1 and e < 1, we have ip^J_^^ ^ * ij;^, hence d^ip'll'^^ < 0. It 
follows that 

'' ifl" ^fo.+)^~dxd^>r,{e,M), 

for M > i? -I- 1, e < 1. At the limit M +oo, e 0, we obtain (gS]). The 
proof of (pS)) is similar. ■ 
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3.2 Proof of Theorem [7] 

Our aim is to show the following 

Proposition 10 Let uo,vo E L°°{M.), < uo,vq < 1 a.e. and let f, resp 
g, be a generalized solution to li^p - jTTp with datum uq, resp. vq- Let M — 
sup |fc(a;)a(^)|. Then we have, for R > 0, 

2:eR,5e[o,i] 



^ [ [ I -f+g-d^dxdt < [ (mo - vo)+dx. (27) 

^ Jo J{\x\<B}JM J{\x\<B.+MT} 

Remark: In case / = Xu, g = Xvj we have J^—f^g-d^ = (u — v)~^ , hence 
(P7)) gives uniqueness of the solution to (|TE| - (fT7|) (more precisely, it gives the 
i^-contraction with averaging in time and the comparison result uq < vq a.e. 
^ u < V a.e.). 

Remark: To obtain the second part of Theorem [71 we apply (^7)) with g — f 
to obtain ^ 

[ [ [ -f+f-d^dxdt < 0. (28) 

^0 J{\x\<R}Jm 

Since < / < xii we have /+ > a.e. and /- < a.e. We deduce from 
(1^51) that /+/- = a.e. Let vt^x denote the Young measure —d^f+: we have 
d^f- — d^f — Sq = and, by examination of the values at ^ = ±oo of /±, 

for a.e. {t, x) £ Q, 

f+{t,x,^) = iyt,x{(„+oo), f^{t,x,^) = -i^tA-oo^O- 

But then, the relation /+/_ — implies that i't,x is a Dirac mass at, say, u{t, x). 
By measurability of v, u is measurable and f — Xu- 

Proof of Proposition llOt Since /+ and g- satisfy 
f+{dtip + k{x)a{S)dx^)dS^dtdx 
+ 11 h,+ij{0,x,Od^dx-{kL-kR)+ I ( a{OHt,0,Od^dt 



d^ipdm+{t,x,£,) (29) 



and 



Q 



g-{dt^ + k{x)a{£,)dxip)d^dtdx 

+ [ [ go,-ij{0,x,Od^dx-{kL-kB)- [ [ a{0i>it,0,0d^dt 
JrJr Jo Jr 

d^i^dp_{t,x,0 (30) 

for all G C^([0,T[xMx]R) (here go,- =sgn_(?;o — ^) andp_ is a non- negative 
measure on [0,T] x R x R supported in [0,r] x R x [0, +oo[), it is possible to 
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obtain an estimate for —f+g^ by setting ip = —g-f in (j29p and "0 = in 
POI) {if being a given test function) and adding the result. This requires first, 
however, a step of regularization. 

Step 1. Regularization. Let Pa,e,s denote the approximation of the unit on 
R'^ given by 

Pa,e,s{t, X, = Pa{t)pe{x)ps{0, 2;, ^ G R^, 

where pe is defined in ((20|) . Let ^ e C^([0, r[xM x R) be compactly supported 
in ]0,T[xR \ {0} x R. Use ^ * Pa,e,5 as a test function in (P^)) and Fubini's 
theorem to obtain 

\ + 



+ / / fo,+i^*Pa,e,siO,x,^)d^dx~{kL~kji)+ J J a{^)'ilj * pa^e,s{t,0,Od^dt 

d^^pdm'^^'\t, X, ^) + Ra,e,s{ip), 



where /"'^''^ /+ * /5Q,e,5, to+^''^' := m+ * pQ^e,^ and 
Ra,E,sW = / / /+[fc(a;)a(C)(c^xV') * Pa,e,5 ~ {k{x)a{£,)dxip) * pa,E,s]d^dtdx. 



Here we have denoted p{t,x,£,) = p(— t, — x, — ^). Also observe that, implicitly, 
we have extended /+ by outside [0, T] since, e.g. 



f+{t)i; * p^{t)dt ^ I / f+{t)^j{s)pa{t-s)dsdt 

i^is) / f+{t)pa{s - t)dtds. 
Jo 

Since ■0 is supported in ]0, T[x]R \ {0} x R, we have, for a, e small enough, 

/o,+V' * (0, X, ^)(i^(i2; =0, 



T 

' a(C)V'*Pa,e,5(<,0,e)def^i =0, 



and 



-Ra,e,<5(0) = f+k{^)[a{i){dxi}) * Pa,e,<5 - {a{S)dx1p) * Pa,e,s]d^dtdx. 

JQ JR 



We deduce 



Q JR 



f+'"\dt^ + k{x)a{i)d^il:)didtdx 



d^Hnil'''\t,x,0 + Rc.,e,5m- (31) 
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A similar work on g_ gives 



= / [ d^^dp^J-''^{t,x,0+Qfi,^A^), (32) 
JQ Jm 

where 

Ql3,u,a{i^) ^ 9-k{x)[a{£_){dxilj) * Pi3,u,cr - {a{£^)dx'ip) * pp^^^d^dtdx. 

JQ JS. 

Step 2. Equation for ^ f'^^'^'^ gPj^^<^ _ Let if e C~([0,T[xR) be non-negative 
and compactly supported in ]0,T[xR\ {0}. Notice that ip does not depend on 
i- Set V = -fg-"'" in (ED), = -<y5/+ in dSS])- Since 

fdt{ipg) + gdtiff) = /^atv? + dti^fg), 

we obtain by addition of the resulting equations 

-f+''^g^-'''''{dtf + k{x)a{Od^(p)d^dtdx 

= - ! V I d^f^''''dp^_'-''''{t,x,0+d^g^_'-'''dmr'\t,x,0 

JQ JR 

Notice that the term 

ip [ d^fr''dp^_^''^-{t,x,£,)+d^gt^'''''dmY\t,x,0 
Q Jm 

is well defined since the intersection of the supports of the functions /"'^'* and 
pf^j^^'^ (^resp. f^''^''^ and m"'^''') is compact. Actually, this term is non-negative 
since p^'^''^ ,m°l''^'^ > and d^f"'^'^ ,d^g^'''''^ < 0. We thus have 

'f+'^'^g-'^''^{dt'P + k{x)a{(,)dx'p)d^dtdx 

> Ra^eA-'fa-"'") + Qp,uA-ff+^'^)- (33) 
It is easily checked that 

hence 

lim Ra.eA-^g-"'") + Qli,v,a{-ff+^'^) = 0" 
0,0-— >0 

At the limit (5, ct ^> in ((33|) , we conclude that 

-f+'g-''{dtf + k{x)a{C)dxip)d£,dtdx > 0. (34) 
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step 3. Traces. Suppose that < k^. We then pass to the hmit £,a — ?> m 
(IMl) to obtam 

f+9-'''{dt^ + k{x)a{^)d^ip)d^dtdx > 0. (35) 

Note that in the opposite case fc^ > kn, and with our method of proof, we would 
first pass to the hmit on /3, v. Let us now remove the hypothesis that Lp vanishes 
at < = 0: suppose that -0 G C^([0,r[xM) is non-negative and supported in 
[0,r[xM\ {0} and apply ^ to ip{t,x) = VX^, a;)Lj^„ (i). We have 

' 1+9-" ^r]„{t){dtip + k{x)a{£_)dx^j)d£,dtdx 

~f+9-'''^Kt, ^W-n„ {t)d^dtdx > 0. (36) 

By (|25|) applied with ip^{x,£_) — g^'''{0,x,^)'4;{0,x), we obtain 



lim / / f+g^j''{0,x,Oil;{0,x)uj' {t)d^dtdx 

> I I /o.+/''^(0,x,OV'(0,x)dfdx. 

JR JR 

Now f+{t, X, £,)9-''{t, X, £,)i^{t, x) has a compact support, say in [0, T] x [—R, R] x 
[— i?, i?], thus ip-{t,x,£,) = g'^'''{t,x,^)'ip{t,x) is uniformly continuous on this 
compact support. Therefore for /i > 0, there exists 7 > such that \ip^ {t, x, ^) — 
(/7_(0, a;, ^)| < for any < i < 7 and any x,£^ £ [—R, R], and then for large 
n, we have ?7„ < 7 and 

/+(i,x,0 [g^J^it^x^Oi^it^x) - g^j''{0,x,0^iO,x)) uj'^Jt)d^dtdx 

8i?2 



< P j P^At) dt ^ iJi. (37) 

Thus we obtain, at the limit n — >■ +00 in p6p . 

'f+g'^'''{dtip + k{x)a{C)dxi')dS,dtdx 

-h,+g^j\G,x,0^{Q,x)didx>G. 



'Q 



The next step is then to remove the hypothesis that ^ vanishes at a; = by 
setting ilj{t,x) — 9{t,x)uj,-i^{x) where 9 G C^([0,r[xR) is a non-negative test- 
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'Q 



function. We have 

'f+g^j''uj^^{x){dte + k{x)a{C)d^e)d£,dtdx 

-f+g'-'^Oit, x)fc(x)a(Ow;^ {x)d^dtdx 

+ 11 -/o,+/^"(0, X, 0^(0, a;)w„„ {x)didx > 0. 
Jm JR 

By m with 0_ (t, = g^J^it, 0, 0), 

hm / / f+k{x)a{£_)iu' {x)g^j'' {t, 0,^)9 {t,0)d£.dtdx 

> -(fci - / / alO/'^^Ci, 0, OOit, 0)d^dt, 

Jo Jm 

and by an argument similar to (1371) . the hmit as [n — > +oo] of the term 

/+fc(x)a(0^;„(x) (^g^_'''{t,x,Od{t,x) - g^j''{t,0,O0{t,0)) d^dtdx 
is zero. We have therefore 

-f+g^lj^idte + k{x)a{^)d^0)d^dtdx 

+ {kL~kR)+ f a{Og^J-'{t,0,O0{t,0)dCdt 

'h+g^j''iO,x,OOiO,x)d^dx >0. 



Since (fc^ — kfj)^ = 0, we have actually 

'f+g^J^idtO + k{x)a{^)d^9)d^dtdx 

'h+g^-''iO,x,OOiO,x)d^dx >0. 



Take /3 = ?7„ where (?7„) is given in Prop. [51 At the limit — > first, then 
n — +00, we obtain 

-f+g-{dtO + k{x)a{C)d^9)d^dtdx 

+ limsup / I -fo,+gl"{Q,x,S)e{Q,x)d^dx>Q. (38) 

Tl-> + 00 jR jR 



Observe that 



f-{0,x,O=f g-{t,x,OPvAt)dt 

T 

g-{t,x,Ouj',Jt)dt. 
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By (transposed to g- tested against a function ip+), we have 

lim / f -fo.+gl"{0,x,00{0,x)d^dx< f f -fo,+go.-0{0,x)d^dx. 

"^+°°JRJR JRJR 

Since 

/ -k.+90-dS, = / -sgn_^(wo - C)sgn_(wo - Od^ = ("o - vo)+ , 
we obtain by (1551) . 

-f+g-idtO + k{x)a{£.)d^e)d^dtdx + [ {uq - wo)+6'(0, a;)dx > 0. (39) 



'Q 

It is then classical to conclude to ^7^: let M > 0,R> MT, let 77 > and let r 
be a non-negative, non-increasing function such that r = 1 on [0, i?], r = on 
[R + rj, +oo[. Set e{t, x) = ^r{\x\ + Mt) in dSH]) to obtain 

^ I I -f+9-r{\x\ + Mt)d^dtdx < [ {uq - vo^dx + J, 

where the remainder term is 

J = / / -f+g-^—^r'i\x\+Mt)iM + k{x)aiOsgiiix))d^dxdt. 
Jq Jm ^ 

By definition of M, J < and since r{\x\+Mt) = 1 for |a;| < R~MT, <t<T, 
we obtain 

Tf; [ [ I -f+g-d£_dxdt < f {uq - vo^dx. 

^ Jo J\x\<R-MT Jm J{\x\<R+rj} 

Replacing Rhy R + MT, and letting 77 -J> gives (gT]). ■ 

4 Convergence of the BGK approximation 

Theorem 11 Let uq G ni°°(M), < uo < 1 a.e. When e ^ 0, the solution 
to the ([ip with initial datum /o — Xuo converges in LP{Q xlS.^), 1 < p < +00 
to Xu, where u is the unique solution to \16\) - ^V?\j . 

Proof: For f (^L'^iRi), set 

™/(0= f {Xu-m)dC, n= f fiOdt 

J —CO Jm 

It is easy to check that mj > if < sgn(^)/(^) < 1 for a.e. ^ (c/. (29) in 
jBre83l ) . In our context, we have < < Xi: hence := -m/e > 0. Viewed 
as a measure, is supported in [0,T] x Mj, x [0, 1]. Integration with respect 
to ^ in dD) gives 

m%0=dt (jy^OdA +dx lk{x) J^iOriOdA 
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in D'(]0,T[xR^). Summing over {t,x) € [0,T] x [xi,X2], C e]0, 1[, wc get the 
estimate 

m^([0,T]x[xi,X2]x[0,l])= / ' / {l-Oir{T,x,O-r{0,x,O)d^dx 



XI Jo 
T rl 



X2 



(40) 



Since f^it) G ^"^(IRk x M^), there exists sequences (a;") I — oo and (x^) t +00 
such that the last term of the right hand-side in PO)) tends to when n — j- +00. 
Since, besides, f^>0 and 

/ \ r{T,x,C)didx < I Xuod^dx ^ \\uq\\l^-k), 
Jr Jr Js. Jr 

we obtain the uniform estimate 

m%[0,T] xMx [0,1]) < ||uo||li(r). (41) 

We also have 

0<,f<Xi, -dJlit,x,0 = ,^tM) + Oie) (42) 

where I't xiO '■— ^u^it,x){0 the identity is satisfied in V {]0,T[xM.x x M^). 
Indeed, by 

r = Xu^ + eidtf + d^{k{x)a{Or)) = Xu^ + 0{e), 

hence 

-d^r+ = -d^f + SoiO = -d^u^ + SoiO + 0{e) = <5„.(0 + 0{e). 

Notice that, for a.e. {t,x), ^ is supported in the fixed compact subset [0,1] 
of M^. We deduce from (|4T|) - (|42|) that, up to a subsequence, there exists a 
non-negative measure m on M"^ supported in [0,T] x Rj, x [0,1], a function 
/ G L°°(]0,T[;Li(M, x R^)) such that < / < Xi, -d^f+{t,x,0 = i^tAO 
where z/ is a Young measure Q — >■ R^ and such that m'^ m weakly in the sense 
of measures (i.e. (m^ — m,ip) — t- for every continuous compactly supported if 
on R3) and / in L°°{Q x R^) weak-star. Besides, since G C([0, T];Ll^^) 

satisfies /^(O) = /o and the BGK equation 

+d^{k{x)a{C)n=d^m', 

it satisfies the weak formulation: for all ^ G Cj?°([0,r[xR x R), 



f%dt^ + k{x)a{^)d^^p)d^dtdx + / / foip{0,x,C)d^dx 

d^ipdrrf{t, x, ^). 
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In particular, we have 

/ / fUdti^ + Hx)a{Od^ij)d^dtdx + f f fo,±HO,x,Od^dx 
JQ Jm Jr Js. 

= / sgn-^{£,)k{x)a{^)dxipd^dtdx + / / d^ipdm^ {t, x, ^) 

Jq jr Jq Jr 

={kR-kL) [ [ sgn^{Oa{O4'{t,0,Od^dt+ [ [ d^^dm'{t,x,0 
Jo JR JQ Jr 

HkL-kn)^ I I aiOHt,0,Od^dt+ [ [ d^^dm'±it,x,0, (43) 

Jo JR JQJR 

where 

- r [ «(e)[(fcL - fcfl)sgn^(e) + (fcL - kn)^m,0,Odidt. (44) 
Jo JR 

More precisely, we set 

m% = m' + a(C)[(fcL - fcH)+sgn+(C) - (fc^ - kii)-8gn_iO]dC5ix - 0), 

and 

ml =m'+ f a{C)[{kL - fc,j)+sgn+(C) - (fcL - fcfl)-sgn_ (C)]dC<5(x = 0). 

J — oo 

Notice that in both cases, and since A{^) > for any ^, we have added a non- 
negative quantity to m^. At the limit e — ;> we thus obtain m^. m± where 
m± is a non-negative measure. Examination of the support of mj_ shows that 
rn+, resp. rri- is supported in [0, T] x x] — oo, 1], resp. [0, T] x Mj. x [0, -|-oo[. 
At the limit e — 0, we thus obtain the kinetic formulation ((T9t . We conclude 
that / is a generalized solution to p^ - pT)) . By Theorem [3 f — Xu where 
u e L°°{Q) is solution to ([TC|) -([T7 | . By uniqueness, the whole sequence (/^) 
converges (in weak-star) to Xu- Actually the convergence is strong since 

/ f \r - Xufd^dtdx ^ f f Iff -2rxu + Xud^dtdx 

JQ JR JQ JR 

< f f r- ^fXu + Xud^dtdx. (45) 

JQ JR 

We have used the fact that < /"^ < 1. The right-hand side of ps)) tends to 
when e — J> since 1, Xu G L°° can be taken as test functions. Hence Xu 
in L'^{Q X M). The convergence in LP{Q x R), 1 < p < +oo follows from the 
uniform bound on in Pi L°°{Q x R). ■ 

Remark: it is possible to relax the assumption that the initial datum for ([T]) 
is at equilibrium and independent on e in Theorem 1111 Indeed, the conclusion 
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of Theorem [TT] remains valid under the hypothesis that the initial datum /q for 
(HI) satisfies 

0</o^<Xi, /o"^/o, uoix):= f foix,Odt (46) 

JR 

where /q ^ /o in (|46)) denotes weak convergence in L^{Rx x R^). Indeed, 
the proof of Theorem [TT] remains unchanged under the following modification: 
passing to the limit in we obtain that / is a generalized solution to 
with an initial datum /o that is not necessary at equilibrium. However, we have 
(c/. (29) in [Bre83]) 

/o - sgn^(^) = sgn±(Mo - C) - 5?™±, 

where m^|_ (resp. mP_) is a non- negative measure supported in [0, T] xRx] — cxd, 1] 
{resp. [0,T] x R x [0,+oo[). Consequently, up to a modification of the kinetic 
measure to±, we obtain that / is indeed a generalized solution to (ITB)) - (IT71) . The 
rest of the proof is similar. 
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